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Abstract 


Lunar and Mars return conditions are examined using the LAURA flowfield code and 
the LORAN radiation code to assess the effect of radiative coupling on axisymmetric ther- 
mochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special 
treatment required to couple radiation in a shock-capturing method is discussed. Results 
indicate that while coupling effects are generally the same as occur in equilibrium flows, 
under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow 
and thus alter relaxation processes. Coupling effects are found to be small for all cases 
considered, except for a five meter diameter aerobrake returning from Mars at 13.6 km/sec. 
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Radiative flux relaxation parameter 


dg 

N Number density, cm -3 

q Heat flux, W/cm 2 

s Path length, cm 

T Temperature, K 

Vco Freestream velocity, km/sec 

z Stagnation line coordinate, cm 

rj Normal coordinate, cm 

Poo Freestream density, kg/m 3 

a Cross-section, cm 2 

r Relaxation time, sec 


Sub- and Superscripts 


e electron 

MW Millikan and White 

n iteration level 

P Park 

R radiative 

s Sound 

t translational 

v vibrational 


Introduction 

A blunt-nosed vehicle traveling through a gaseous medium at hypersonic speed may 
experience significant radiative heating. This heating results from radiative energy emitted 
in the shock layer without directional preference. That portion of energy which is radiated 
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away from the vehicle and escapes the shock layer causes the phenomenon called radiative 
cooling. So far, flows with significant radiative cooling have been simulated only by fixed- 
shock schemes including the viscous-shock-layer (VSL) scheme 1-4 or inviscid equilibrium 
schemes 5 . This approximation simplifies the calculation of the coupled, radiating flowheld 
but may not be valid for flow situations characteristic of aerobraking reentry into the Earth’s 
atmosphere upon return from the Moon or Mars. 

The objective of the present work is to examine radiative coupling effects in shock- 
capturing, thermochemical nonequilibrium flowheld calculations. A number of freestream 
conditions, representative of Lunar and Mars return proposals, are studied. The Langley 
Optimized RAdiative Nonequilibrium (LORAN) code 6,7 , is coupled to the Langley Aerother- 
modynamic Upwind Relaxation Algorithm (LAURA) howheld solver for this purpose. LO- 
RAN includes options for both tangent slab and an approximate 3-D radiation transport 
calculation. Results presented in this paper utilize the tangent slab approximation. LAURA 
is a shock-capturing Navier Stokes solver which includes 11-species chemical nonequilibrium 
and the two-temperature thermal nonequilibrium model of Park 8 . Coupling of the LAURA 
and LORAN models is achieved through an iterative procedure, in which the radiation prop- 
erties are updated every 200 to 5000 howheld iterations, depending on the how conditions. 
Special treatment of the shock and shock precursor regions is found to be necessary and is 
described here. Additional numerical considerations are also necessary to establish a con- 
vergent, efficient scheme. 

The coupled LAURA/LORAN code has been applied to a range of axisymmetric how- 
helds representative of missions of current interest. The effects of radiation coupling in a 
shock-capturing scheme for thermochemical nonequilibrium howhelds with signihcant ion- 
ization are found to differ from those observed in equilibrium or near-equilibrium howhelds. 
In particular, radiative cooling can alter the ionization chemistry in the howheld, resulting 
in an increase in translational temperature and shock standoff distance. 
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Radiation Coupling Theory 


Radiation is inherently coupled to the state of the gas from which it emanates. It 
results from radiative transitions in the internal energy states of the atoms, molecules and 
ions present in the gas. The radiation emitted by such transitions travels through the gas 
at the speed of light until it is reabsorbed or exits the flowheld. Radiation thus provides an 
additional energy transport mode. Radiation can also affect the chemical state of the gas 
through phenomena such as photodissociation, photoionization, and radiative recombination. 

Energy transport through the gas by radiation is modeled in the Navier- Stokes equa- 
tions with a term which appears in the energy equation: the divergence of the radiative 
flux, V • q-R (denoted Q ra d by Gnoffo 9 ). This divergence term is computed in the LORAN 
radiation-transport model. When the tangent slab approximation is used, V ' = dq^/dq. 

The radiation effects on chemical kinetics can be included by adding photo-reactions to 
the chemistry model of the gas and accounting for their effect on the radiation properties. 
The radiation effects on chemical kinetics are not included in the current LAURA/LORAN 
code. This approximation worsens as the energy in the flow increases, but has not been 
quantitatively assessed. 

The transport of energy due to radiation is usually neglected in the aerothermody- 
namic analysis of hypersonic flowhelds. For flowhelds where this assumption is not valid, 
the effects of radiation can be determined by comparing solutions from calculations with 
and without radiative transport. In general, the shock layer becomes non- adiabatic when 
radiation transport is included. The enthalpy of a parcel of fluid decreases as it traverses the 
layer from the shock to the body due to radiative cooling. For equilibrium gas calculations, 
the resulting temperature profiles are no longer constant but decrease monotonically across 
the layer. Lower average temperatures result in higher average densities and thus decreased 
shock standoff distance. The combination of lower temperatures and a thinner shock layer 
reduces the radiative heat load to the wall. 

In a thermochemical nonequilibrium gas calculation, multiple temperatures are defined 
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to describe the partitioning of energy across the internal modes. The number of temperatures 
required to adequately describe the flow depends on the conditions. For example, rotational, 
vibrational, and electronic energy modes for the mixture or for each species may be used. 
The various states of each energy mode are assumed to be in equilibrium at that mode’s 
characteristic temperature. Each mode is described by its own energy conservation equation 
which includes terms coupling it to the other energy modes. Usually, separate internal modes 
for each species are not defined since the energy modes are strongly coupled between species. 

To correctly model the radiative coupling in the nonequilibrium case only that portion 
of V • qR associated with transitions of each particular energy mode should be included in 
the energy equation for that mode. Since purely rotational radiative transitions occur at 
low infrared energies — which are generally neglected in gas radiation — it may be assumed 
that all radiative transitions of interest in air involve a transition between vibrational or 
electronic states. LAURA employs the two-temperature model 8 , requiring a total energy 
equation (T t ) and a combined vibrational - electron - electronic energy equation (T v ). Thus 
the complete V • qn appears in both energy equations. 

Numerical Algorithm 

The LAURA algorithm for solution of the Navier Stokes equations is summarized in 
Ref. 10. The divergence of the radiative flux can be treated as a source term in the two 
energy equations. In general, source terms are handled implicitly in the LAURA algorithm. 
This treatment improves the scheme’s stability and robustness. Since radiation is sensitive 
to flowheld properties, it would be advantageous to treat the added radiation source term 
implicitly also. Park and Milos 11 , however, have found this implicit treatment to be pro- 
hibitively expensive (in computer time) for anything except 1-D flows. Furthermore, it is 
not clear how to obtain the necessary terms for the Jacobian in the current nonequilibrium 
treatment. The radiation term is therefore treated as an explicit source term in the LAURA 
algorithm. 
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The convergence of the flowheld in LAURA is characterized by transient waves in the 
gas properties. When the sensitive radiation calculation is performed on a flowheld still 
containing transients, the radiation term can change dramatically. This can introduce addi- 
tional fluctuations and lead to instability. In the present study, therefore, under-relaxation 
is applied to the explicit radiation source term. That is, the latest radiation solution is 
combined with the previous solution to compute an updated radiation term: 

(v • <m)" ,+1 = d q {v- m ) 71 + (i - d q ) (v • qR ) n+1 (i) 

The under-relaxation approach of Eq. 1 reduces the fluctuations due to coupling and increases 
the stability of the combined algorithm. A relaxation factor, d g , of 0.5 is used for the results 
presented. 

Because the how properties change little in a single howheld iteration, and because 
it is advantageous to allow howheld transients to decay, it is not necessary to update the 
radiative hux divergence term every time the howheld is updated. In fact, the total cost 
of the combined computation can be reduced several orders of magnitude by updating the 
radiation term a minimum number of times, as is successfully done in equilibrium shock- 
fitting methods. Not updating it frequently enough, however, can decrease the convergence 
rate, produce oscillatory non- convergent behavior, or destabilize the solution. For example, 
if the howheld changes dramatically between radiation term updates, a large change in the 
radiation will result and can introduce severe waves into the howheld. In addition, if how 
structures change position signihcantly between updates of the radiation, the radiation term 
will temporarily not be “aligned” with the associated how structure. Experience with cases 
characterized by a moderate degree of nonequilibrium suggests that updating the radiation 
calculation once for every 500 howheld updates provides reasonable convergence. For cases 
with more severe nonequilibrium in which coupling has a large effect on the standoff dis- 
tance, more frequent updates are required (every 200 iterations). As a howheld approaches 
equilibrium, on the other hand, the current LAURA algorithm converges more slowly. For 
near-equilibrium cases, the radiation is updated only every 1000 to 5000 howheld iterations. 
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Obtaining coupled solutions may be difficult when extreme gradients occur in the shock 
wave, as is the case for near-equilibrium flows. The computational mesh used in a flowheld 
calculation seldom provides adequate resolution in the shock region to produce smooth tem- 
perature profiles in such flows. This lack of resolution at the shock does not significantly 
degrade the flowheld solution. Since the radiation calculation is very sensitive to tempera- 
ture, however, it can cause large changes in the profile of dq^/dq which may lead to inaccurate 
results or even divergent solutions. Increasing the number of grid points can alleviate this 
problem but is not always practical because of storage limitations. Grid adaption, on the 
other hand, can improve the temperature profiles without increased storage. To provide this 
ability, a grid adaption algorithm has been added to the LAURA code. In this study it is 
used to cluster grid points in the vicinity of the maximum temperature in the shock. 

The low-temperature, or precursor, region ahead of the shock in a shock-capturing 
scheme also requires special treatment. The LORAN radiation model as currently imple- 
mented does not include accurate radiation properties at low temperatures. It is generally 
assumed that the precursor has a negligible effect on the converged flowheld, since it can ab- 
sorb only a small fraction of the radiative energy passing through it. Recent results 12 confirm 
the validity of this approximation. Therefore, dq^/dq is set to zero - the value appropriate 
for a non-participating medium - when either temperature (vibrational or translational) is 
below a certain minimum. This is done within the LAURA howheld calculation every itera- 
tion. Even if accurate precursor (low - temperature) radiation properties are available, some 
logic is still required so that any decrease in shock standoff distance which might occur be- 
tween radiation updates will not leave a large radiation term in what is now the freestream. 
The value used for this minimum temperature is 500 K for the cases in this paper. 

Physical Models in LAURA 

The complete set of physical models used in the LAURA code is described in Ref. 9. Op- 
tions exist for some of these models and a number of parametric studies have been performed 
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to assess the differences between them 13 16 . For the results presented here, LAURA utilizes 


the Park-91 chemical kinetics model 17 (except for the Mars capsule return case which uses 


the Park-87 model 8 ). The dissociation rate controlling temperature is modeled by \jT t T v 
for the LTV and Mars aerobrake cases, and by TfT^ 3 for the capsule return cases. The 
vibrational - translational relaxation time is modeled as 


^ MW , P 

Tv = T„ + T 


with 

rf = (oyUlV ) -1 

where the effective cross section for vibrational relaxation is 


( 2 ) 


( 3 ) 


a v = 10“ 17 (50, 000/Tt) 2 cm 2 


( 4 ) 


and t„ p is calculated with the total number density. Relative to results presented earlier 13 , 
this causes faster relaxation of the vibrational energy mode (as was suggested in Ref. 13 
by comparison to the FIRE II results). The vibrational energy removed in dissociation is 
assumed to be the average vibrational energy, while in ionization it is assumed to be the 
energy to ionize from the ground state. 

The wall boundary condition for the LTV and Mars aerobrake cases is Stewart’s finite 
wall catalysis model 18 . The capsule return cases employ the “super-catalytic” boundary 
condition in which freestream species concentrations are enforced at the wall. 

While these models have been validated to the extent possible, using FIRE II flight 
data and ground-based experimental results, they nevertheless retain significant uncertainties 
as does the radiation model. Their application to the Mars return conditions represents 
an extrapolation from any known data. The results presented below should be viewed as 
presenting trends rather than absolute heating levels. 


Results 


The effects of radiation coupling on thermochemical nonequilibrium flows are illustrated 
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by comparing uncoupled and coupled solutions from the LAURA/LORAN code. An uncou- 
pled wall radiative flux is computed by applying the LORAN code to the uncoupled LAURA 
output. Flow conditions representative of aerobraking reentry of axisymmetric sphere cone 
bodies into the Earth’s atmosphere upon return from the Moon and Mars are examined. 

Lunar Return Conditions 

Two lunar return conditions are studied. The first is a small body representative of a 
capsule return. The peak convective heating point on this trajectory at 9.871 km/sec and 
75 km altitude ( p ^ = 4.623 X 10 -5 kg/m 3 ) is examined. A 42 by 64 cell grid defines the 
forebody of this 3 m nose radius, 70° sphere cone capsule. The second case is representative 
of a large lunar transfer vehicle (LTV). The peak convective heating point at 9.797 km/sec 
and 72 km altitude ( p ^ = 5.987 X 10 -5 kg/m 3 ) is computed. A 30 by 64 cell grid defines 
the forebody of this spherical aerobrake which is 15.24 m in diameter, with a 13.58 m nose 
radius and a 0.3048 m radius shoulder. 

Lunar Capsule Return 

The stagnation line temperature profiles from the uncoupled solution for the capsule 
return case are shown by the solid lines in Fig. 1. Thermal nonequilibrium effects are 
restricted to about 20 percent of the shock layer in this near-equilibrium case. The dashed 
lines in the figure are the temperature profiles from the coupled solution. For this near- 
equilibrium case, the behavior of the coupled solution is the same as has been observed with 
equilibrium flowhelds: the radiation cools the flow, leading to lower temperatures (which do 
not level off to an equilibrium value as in the uncoupled solution), higher density, and thus 
a smaller shock standoff distance. The result is a temperature which is 200 to 250 K lower 
through much of the shock layer. There is no detectable effect on the vibrational relaxation 
distance: the temperatures relax to thermal equilibrium in about 3 cm whether coupling is 
included or not. 

The effect of coupling on the surface heating contributions is shown in Fig. 2. Radiative 
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cooling results in lower surface radiative heating all along the body. At the stagnation point 
the coupled result is 7 percent lower. It is 20 percent lower on the body flank, where the 
volume of radiating gas in the shock layer is reduced the most for this sphere cone geometry. 
The effect on convective heating is small. It decreases a maximum of three percent at the 
stagnation point. This insensitivity is expected because of the competing effects of lower 
shock layer temperatures, boundary layer absorption of radiation, and thinning of the shock 
layer. 

For the near-equilibrium capsule return case, starting from a converged uncoupled solu- 
tion, about 22000 flowheld iterations with eight radiation updates were required to achieve 
the converged coupled solution. This required about eight hours of Cray-2 CPU time. The 
reason for this large computer time is the very near-equilibrium character of this case. In 
order to provide a smooth solution for the wall radiative heat flux, the temperature gradients 
in the shock must be resolved to a much greater accuracy than is required for a reasonable 
flowheld solution. For this case, the necessary resolution was achieved through a trial and 
error succession of grid adaptions, convergence runs, and radiation updates. An optimal 
selection of the adaption and radiation update intervals would have required less CPU time. 

Lunar Transfer Vehicle 

The stagnation line temperature profiles for the LTV are shown in Fig. 3. This case is 
very near equilibrium, with thermal relaxation occurring in the outer 7 percent of the shock 
layer. The effect of radiative coupling is small. It reduces the shock standoff distance by 
about 4 percent and produces continually decreasing temperatures through the shock layer. 
The temperatures are 300 to 400 K lower over a substantial part of the shock layer. These 
trends are again the same as for equilibrium. 

The effect on the surface heating for the LTV is shown in Fig. 4. For this spherical 
forebody, the radiative flux decreases monotonically away from the stagnation point. The 
radiative heat flux is reduced 15 to 20 percent by coupling, while the convective heating 
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decreases 7 to 9 percent. The larger radiative flux for this large body results in increased 
radiative cooling effects. 

For this near-equilibrium case, starting from a converged uncoupled solution, about 
22000 flowfield iterations were again required to obtain the coupled solution. Four radiation 
updates were performed. Fewer radiation updates were required for this case because the 
shock moved very little with coupling. The calculations required just over three hours of 
Cray Y-MP CPU time for this smaller grid. 

Mars Return Conditions 

Four Mars return conditions are considered. The first is a small body representing a 
capsule return. The flight condition is 12 km/sec and 80 km altitude ( p ^ = 1.999 X 10 -5 
kg/m 3 ) with a fine 30 by 128 cell grid covering only the nose region of a 1.08 m nose radius, 60° 
sphere cone. This case provides a detailed look at the stagnation region of this flowfield. The 
remaining three cases involve a geometry and conditions typical of an aerobraking reentry 
into the Earth’s atmosphere 19 . The geometry is a scaled down version of the LTV which is 5 
m in diameter with a 4.46 m nose radius and a 0.1 m radius shoulder. The grid is 30 by 64 
cells. The three freestream conditions studied here span the entrance trajectory from near 
the continuum limit down to the predicted peak convective heating point. 

Mars Capsule Return 

The stagnation line temperature profiles from the uncoupled solution for the Mars 
capsule case are shown by the solid lines in Fig. 5. Half the shock layer is in thermal 
nonequilibrium. Temperature profiles for the coupled solution are also shown in Fig. 5 by 
the dashed lines. The effects of coupling are small but still affect the relaxation process in 
this nonequilibrium flow. While the vibrational temperature is again reduced by radiative 
cooling, the translational temperature increases in the nonequilibrium region. This unex- 
pected behavior is attributed to thermal nonequilibrium effects in significantly ionized flows 
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(25 percent ionization here vs 4 percent in the lunar cases). Examination of the electron 
number density profiles (Fig. 6) shows that the electrons are freed in two discernible surges: 
one immediately behind the shock due to the associative ionization reactions and one deeper 
in the shock layer due to electron impact ionization. The associative ionization reactions 
provide the free electrons necessary to initiate electron impact ionization. Since the vibra- 
tional temperature in this region is reduced as much as 800 K by coupling, both types of 
ionization are slowed. As a result the electron number densities are up to 75 percent different 
in the region of maximum coupling effects, and the maximum number density of electrons 
in the shock layer is reduced by 12 percent. The result is that the translational temperature 
remains higher (less energy is converted by ionization). The density in the shock layer there- 
fore remains lower in the nonequilibrium region. With the compensating effect of higher 
density near the wall, the net effect is no change in the shock stand-off distance. 

Radiation cooling has a small effect on the wall heat fluxes for this condition as shown 
in Fig. 7. The radiative flux is 5 to 9 percent lower in the coupled case, while the convective 
heating is 6 to 10 percent lower. (Recall that only the nose region of this aerobrake flowheld 
is computed.) 

The convergence of this relatively nonequilibrium case was much smoother than for 
the near-equilibrium cases. Starting from a converged uncoupled solution, it required 4250 
LAURA iterations and five radiation updates. This was done in under two hours of Cray 
Y-MP CPU time. A significant savings in this case is that grid adaption is unnecessary to 
resolve this diffuse shock. Once the initial transients from radiation were absorbed in the so- 
lution, radiation was updated every 500 LAURA iterations. Examination of the convergence 
histories suggests that this is a near-optimal interval. 

Mars Aerobrake Return — Point 1 

To examine the behavior of thermochemical nonequilibrium coupled flows, a point early 
in the aerobraking trajectory is studied. The velocity is 13.604 km/sec at 91.5 km altitude 
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( Pco = 1.239 x 10- 6 kg/m 3 ). 

Stagnation line temperature profiles for this case are shown in Fig. 8. The very nonequi- 
librium character of this flow is evident; the temperatures are out of equilibrium through 
most of the shock layer. The coupled profiles (dashed lines) show qualitatively the same 
behavior observed for the Mars capsule case, but with more pronounced effects. Radiation 
cooling decreases the vibrational temperature enough to inhibit ionization reactions in this 
25 percent ionized flow. This results in more energy staying in the translational energy 
mode and a longer relaxation time. (Dissociation of molecules is also inhibited by the lower 
vibrational temperature, but this has less impact on the flow because the number densities 
and energies of formation are lower than for the ionized species.) The density is lowered 
and the shock standoff distance increases. The electron number density profiles are shown 
in Fig. 9. They show that the electron avalanche due to electron impact ionization reactions 
is delayed by about 10 cm inside the shock layer. As a result the maximum electron number 
density in the shock layer is reduced by a factor of three from the uncoupled result and the 
ionization fraction is lowered to 20 percent. 

The effect on surface heating is shown in Fig. 10. Though the radiative heat load is 
small, nonequilibrium coupling effects are strong enough to reduce it by a factor of two or 
more. The convective heating increases 35 to 70 percent with coupling, due to reabsorption 
of radiation in the boundary layer. This is indicated by the higher temperatures near the 
wall for the coupled solution (Fig. 8). Despite the increase in convective heating, the total 
wall heat flux is reduced 18 to 25 percent by radiation coupling effects. Furthermore, while 
the uncoupled solution predicts radiation to be twice as strong as convective heating, the 
coupled solution has convection as the dominant mode of heat input to the wall. 

Convergence of the coupled result for this very nonequilibrium case required more fre- 
quent radiation updates due to the large motion of the shock and the strong coupling effects. 
Updating radiation every 500 to 1000 iterations was found to cause an oscillatory behav- 
ior which converged very slowly if at all. To obtain reasonable convergence rates, it was 
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necessary to update radiation every 200 LAURA iterations. Including the early oscillatory 
problems, convergence of this case took 9000 LAURA iterations, with 19 radiation updates. 
The complete trial and error process took under four hours of Cray Y-MP CPU time. 

Mars Aerobrake Return — Point 2 

A point 30 seconds later in the trajectory is studied to investigate radiative coupling 
effects as equilibrium is approached. Here the velocity is 13.585 km/sec and the free-stream 
density is about an order of magnitude higher at 1.721 xlO -5 kg/m 3 (80.7 km altitude). The 
nonequilibrium region is about a quarter of the shock layer, as shown by the temperature 
profiles in Fig. 11. Coupling reduces the stagnation point shock standoff distance by about 
ten percent. The relaxation zone is not qualitatively affected, however, and remains about 
5 cm long. 

The effect of the cooling on the wall radiative flux is large, resulting in a 25 to 35 percent 
reduction as shown in Fig. 12. The convective heating is unaffected by radiative coupling. 

Converging this coupled solution starting from the uncoupled flowheld required 9600 
LAURA iterations with six radiation updates. For this nonequilibrium case, radiation was 
generally updated every 500 to 1000 LAURA iterations. This solution required just under 
two hours of Cray Y-MP CPU time. 

Mars Aerobrake Return — Point 3 

The peak convective heating point (12.883 km/sec, p ^ = 1.598 X 10 -4 kg/m 3 at 65 km 
altitude) on the aerobraking trajectory is an equilibrium flowheld. The relaxation zone is 
thinner than that shown in Fig. 3. The nonequilibrium LAURA/LORAN code is not well- 
suited to this equilibrium problem. This case can, however, be used to examine the accuracy 
of the nonequilibrium LORAN code at equilibrium conditions. 

Radiation was computed for the converged uncoupled LAURA howheld with nonequi- 
librium LORAN and with equilibrium LORAN. The equilibrium version of LORAN replaces 
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the quasi-steady-state (QSS) 20 excited state populations with Boltzmann distributions. The 
distributions of wall radiative heat flux are compared in Fig. 13. They differ by two to 
eight percent. For an independent comparison, coupled and uncoupled inviscid equilibrium 
stagnation point solutions from the method of Ref. 5 are also shown in this figure. The 
agreement between the three uncoupled solutions at the stagnation point verifies that the 
nonequilibrium LORAN code approaches the equilibrium limit. The coupled inviscid equi- 
librium stagnation point solution indicates the large extent of radiative cooling present in 
this flow. 

Summary of Results 

The results from these six cases can be summarized in a form previously used to estimate 
coupling effects in equilibrium flows. Fig. 14 shows the cooling factor, q CO upied/ ([uncoupled , as a 
function of the uncoupled radiation energy fraction, q U ncoupied/ (PooV^/2). The denominator 
in the latter quantity is a measure of the total energy in the freestream. The broad cross- 
hatched region represents extensive calculations with the equilibrium method of Ref. 5. The 
six nonequilibrium cases calculated in the present work are overlaid on this plot as circles. 
While they do not fall entirely within the equilibrium band, the nonequilibrium cooling 
factors show general agreement with the previous equilibrium results. 

Discussion 

Both LAURA and LORAN are collections of numerous, often complicated, physical 
models. There are questions about the validity of some of these models at the conditions 
examined. In particular, LAURA is a continuum Navier-Stokes model which assumes a 
weakly ionized two-temperature gas. The Mars return cases include degrees of ionization as 
high as 25 percent, which may violate the weakly ionized assumption. Point 1 for the Mars 
aerobrake case is at 91.5 km altitude - near the edge of the continuum region. In LORAN, 
the QSS assumption is open to question in high gradient captured shocks. Furthermore, 
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neglecting the influence of radiation on photodissociation and photoionization may not be 
valid at these extreme conditions. 

The reported computational cost of the coupled calculations is high in part because of 
uncertainty in the optimal frequency for updating the radiation calculation. This optimal 
frequency is shown to be related to the degree of nonequilibrium in the flowheld. Experience 
will lead to insight into selecting this frequency and may reveal an automatic mechanism 
which will result in a more efficient combined LAURA/LORAN code. 

The results presented in this paper use the tangent slab transport approximation. It 
has been noted before 6 that the variation of V ' <?R obtained from the approximate 3-D 
transport solution is smoother than the variation of dq^/dq obtained from the tangent slab 
model. This smoothness would eliminate the need for some of the numerical considerations 
described above but at present represents a significant increase in computational cost. Future 
work will examine coupling with 3-D radiative transport. 

The results suggest that thermochemical nonequilibrium coupled radiation effects are 
important only at very low density (i.e. high altitude), high velocity conditions. The tra- 
jectories of most Earth entry vehicles continue deeper into the atmosphere so that the high 
altitude heating is of minimal importance compared to the heat loads encountered at the 
lower altitudes. If an aerobraking trajectory is proposed which dissipates its kinetic energy 
in the upper, rarefied regions of the atmosphere, significant radiative cooling effects must be 
expected. 


Conclusions 

The LORAN thermochemical nonequilibrium radiation code is coupled with the LAURA 
flowheld code and applied to howhelds characteristic of Earth entry upon return from the 
Moon and Mars. Coupled solutions are obtained using the tangent slab radiation transport 
model. Coupling a radiation calculation to a howheld solution by an iterative procedure 
gives rise to potential instabilities. These are avoided by the use of under-relaxation. In 
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addition, special treatment of the radiation is required in the vicinity of the shock and in 
the shock precursor region which is not necessary in shock-fitting methods. 

The results of this study indicate that radiation coupling effects are small for the lu- 
nar return cases considered. For these near-equilibrium cases, coupling effects are similar 
to those predicted by equilibrium shock-fitting methods. While the combined nonequilib- 
rium LAURA/LORAN code is not an efficient means of computing equilibrium flowhelds, 
it is shown to produce the correct equilibrium limit. For the higher energy Mars return 
cases, both nonequilibrium and coupling effects are stronger. In these cases with significant 
ionization, radiation cooling can alter the ionization chemistry of the flow. This change 
in chemistry can produce increased translational temperatures and result in a larger shock 
standoff distance. These effects are unlike those observed at equilibrium or near-equilibrium 
conditions. Nevertheless, the cooling factors for all these cases show agreement with previous 
equilibrium experience. 
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Figure Captions 


1. Effect of radiation coupling on stagnation streamline temperature profiles for the 
lunar capsule 

2. Effect of radiation coupling on surface heating for the lunar capsule 

3. Effect of radiation coupling on stagnation streamline temperature profiles for the 
lunar transfer vehicle 

4. Effect of radiation coupling on surface heating for the lunar transfer vehicle 

5. Effect of radiation coupling on the stagnation streamline temperature profiles for the 
Mars capsule return condition 

6. Effect of radiation coupling on stagnation streamline electron number density profiles 
for Mars capsule return condition 

7. Effect of radiation coupling on surface heating for the Mars capsule return condition 

8. Effect of radiation coupling on the stagnation streamline temperature profiles for the 
Mars aerobrake return condition - Point 1 

9. Effect of radiation coupling on electron density for the Mars aerobrake return condi- 
tion - Point 1 

10. Effect of radiation coupling on surface heating for the Mars aerobrake return con- 
dition - Point 1 

11. Effect of radiation coupling on the stagnation streamline temperature profiles for 
the Mars aerobrake return condition - Point 2 

12. Effect of radiation coupling on surface heating for the Mars aerobrake return con- 
dition - Point 2 

13. Wall radiative heating predictions for the Mars aerobrake return condition - Point 
3 

14. Cooling factor 
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Figure 1: Effect of radiation coupling on stagnation streamline temperature profiles for the 
lunar capsule. 
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Figure 3: Effect of radiation coupling on stagnation streamline temperature profiles for the 
lunar transfer vehicle. 



Figure 4: Effect of radiation coupling on surface heating for the lunar transfer vehicle. 
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Figure 5: Effect of radiation coupling on the stagnation streamline temperature profiles for 
the Mars capsule return condition 



Figure 6: Effect of radiation coupling on stagnation streamline electron number density 
profiles for Mars capsule return condition. 
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Figure 7: Effect of radiation coupling on surface heating for the Mars capsule return condi- 
tion. 



Figure 8: Effect of radiation coupling on the stagnation streamline temperature profiles for 
the Mars aerobrake return condition - Point 1. 
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Figure 9: Effect of radiation coupling on electron density for the Mars aerobrake return 
condition - Point 1. 
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Figure 10: Effect of radiation coupling on surface heating for the Mars aerobrake return 
condition - Point 1. 
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Figure 11: Effect of radiation coupling on the stagnation streamline temperature profiles for 
the Mars aerobrake return condition - Point 2 



Figure 12: Effect of radiation coupling on surface heating for the Mars aerobrake return 
condition - Point 2. 
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Figure 13: Wall radiative heating predictions for the Mars aerobrake return condition - Point 
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Figure 14: Cooling factor 
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